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Abstract 

Solutions of equations of geodesic deviation in three- and four-dimensional 
spaces are obtained via the inverse scattering transform. It is shown that in 
the case of three-dimensional space solutions of geodesic deviation equations 
are reduced to solutions of the well-known Zakharov-Shabat problem. In four- 
dimensional space system of geodesic deviation equations is associated with a 
3x3 matrix Schrodinger equation, and dependence on parameters is defined 
f"^ ■ by the nonlinear equations of three- wave interaction. 

Keywords: geodesic deviation, mKdV equation, Chandrasekhar metrics, matrix 
Schrodinger equation. 



1 Introduction 

It is "well-known that a m x m matrix Schrodinger equation on — oo < x < oo is 
defined by the following expression £Q : 

O- Lip(x,k) = Xtp(x,k), A = fc 2 , 

where 

X ' L=-(d 2 /dx 2 )I + U{x), 

1 = (<%)> U(x) = (uij(x)) ; i,j = 1, ...,m, 
ip(x,k) = [ipi(x,k),ip2(x, k), . . .,ip m (x, k)] . 

Further, let rf be the components of deviation vector between two infinitesimally 
nearby geodesic lines. Then the components rf satisfy to the Jacoby equation |2] 

jVi{r?V s r}) = -v l R\ km v m rf, (1.1) 

where v l are the components of the tangent vector along a geodesic line 7, is 
the curvature tensor of the metrics 

ds 2 = gijdx l dxK 

In a special system of coordinates, where axis x J is a geodesic line, equation l|l.l|) 
has the following form 

^ + 4^=0- (1.2) 

In the paper |Sj it has been shown that in the case of three-dimensional space 
with the metrics 

ds 2 = dx 2 + A(x, y, z)dy 2 + 2B(x, y, z)dydz + C(x, y, z)dz 2 (1-3) 
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the equations of geodesic deviations 



d 2 n 2 

+ R 2 12lV 2 + R 2 131 if = 0, 
^ + R 3 12lV 2 + R 3 13lV 3 = (1.4) 
may be represented in the form of the 2x2 matrix Schrodinger equation 

+ (~Rl2i)v 2 + (-Rlsi + A 2 )r/ 3 = \W- (1-5) 
On the other hand, it is known that AKNS-system jS] 

^L + iA-0i = q(x,y,z)il>2, 

= r(x,y,z)ipi (1.6) 
can be rewritten in the form of a Schrodinger-like equation 



1 \ I rq q x 

1/ 9x 2 I r x rq 



The comparison of the systems <|1.7[1 and l|1.5|l gives the following conditions on the 
curvature tensor 

\ 2 -R 2 21 =rq, \ 2 ~ R 3 131 = rq, 



Analogously, in the case of 4-dimensional space with a geodesic coordinate sys- 
tem 

ds 2 = di 2 + g ab dx a dx b (1.9) 
the geodesic deviations equation has the form 

■ ' RoioV 1 + R^ov 2 + Rlwv 3 = o, 



dt 2 

d 2 r, 2 

~~dW 

d 2 rj 3 
~dt 2 ~ 



+ RUv 1 + R 2 Q 2oV 2 + R 2 o 3 oV 3 = 0, (1.10) 
+ RmoV 1 + Rl 2Q V 2 + RlsoV 3 = 0. 



In the present paper we consider solutions of the equations i|1.4|) and (|1.10(l ob- 
tained by the inverse scattering transform. Our consideration is realized on the basis 
of a Chandrasekhar metrics |H] (the so-called space-time of a sufficiently general 
structure), which includes as particular cases the static and spherically symmet- 
ric solutions (Schwarzschild and Reissner-Nordstrom metrics), and also stationary 
and axially symmetric solutions (Kerr and Kerr-Newman metrics) and so on. In 
section 2 we introduce a three-dimensional analog of the Chandrasekhar metrics, 
the particular case of which is coincide with the metrics (jl.31) . It is shown that in 
the orthonormal basis, related with this metrics, solutions of the system ljl.4JI are 
reduced to the solutions of the Zakharov-Shabat problem 0. Thus, a dependence 
of the potential u on parameters y and z is described by the modified Korteweg-de 
Vries (mKdV) equations. Different particular cases, in which the vector of geodesic 
deviation r/ is explicitly expressed via the fundamental solutions (Jost functions) of 
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the Zakharov-Shabat problem, are considered at the end of section 2. In section 3 
we introduce a 3 x 3 matrix Schrodinger equation which then is associated with the 
system of type (jl.lOf) . Further, a dependence on parameters is reduced to evolution 
equations of the well-known problem of three- wave interaction, the explicit solutions 
of which was obtained by Zakharov and Manakov in 1973 |10l 1111 ITS] . It is shown 
that in the case of decay instability and reality of potential matrix, the system of 
equations of geodesic deviation (|1.10|l has a wide class of particular solutions. 

2 Three-dimensional space 

2.1 The three-dimensional Chandrasekhar metrics 

Let us consider in the three-dimensional space with a signature (—,—,—) a metrics 
of the following form 



ds 2 - - ]T e 2 ^ A (dx A ) 2 - e 2 ^(dx 3 - <Mcb 



A\2 



(2.1) 



where A = 1,2. ip, ^A and qA are the functions on variables x , x 2 , x 3 . 

The orthonormal basis, related with this metrics, is defined by the following 
covariant and contravariant vectors 

e (1)i = (0, 0, -e' 11 ), e (a)i = (0, -e^, 0), 



e(3)i = (-e^, qie*, q%e^). 
= fee - " 1 . °> e (2) = (<?i^ 2 , e-»\ 0), 



From ((221 and (O it follows that 



(e-* 0, 0). 





-1 








V(a)(b) - 





-1 













-1 



Let 



(2.2) 



(2.3) 



(2.4) 



lj a = e fiA dx A , u 3 = e^(dx 3 ~ ^ q A dx A ) 

A 

be the basis 1-forms. It is easy to see that inverse relations for (|2.4[) have the form 
dx A = e-^ A uj A , dx 3 = e~^u) 3 + ^ e-^aA (2.5) 



Expressing the exterior derivatives of the forms u> 1 via the basis 2-forms 
w'a^ (i ^ j, i,j = 1, 2, 3), we have 



duj A = ^ e^ A ha,b<1x b /\dx A + e^ A ha,3<1x 3 Ada 

B 

= Yj e ~ M " ^ A > B + <1bHa,3)u B Au A + e^M, 3 w 3 Aw A (2.6) 



e^u 3 
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For the brevity of exposition let us introduce a derivative of the function f{x l , x 2 , x 3 ) 
on a coordinate x A (A= 1,2) which we will denote as f :A , 

f:A = f.A + q A f, 3 . (2.7) 
This operation is the differentiation, since it satisfies to a Leibnitz rule 

(fg)-.A = fg-.A + gf-.A- 

Using H2.7[) . we can rewrite the equation l|2.6|l in the form 

&u A = - J2 e-^HA:BU A A" B - e-'V^'W. (2.8) 

B 

In like manner we have 

do, 3 = Y, (V^a + q A ,z)u A h^ - Y e^ A -» B q A :BU A hu B . (2.9) 

A A,B 

Further, the equations 



= du) 1 +ivjf\uj l = (2.10) 

1 



^Rj km u; k Au m = Slj (2.11) 

are called respectively the first and second Cartan structure equations, where the 
Cartan 2-form Qj is 

n{ = duj + wjAwf • 
Owing to absence of torsion (T J = 0) the first Cartan structure equation gives 

du> 3 = - ^u^AuA (2.12) 

A 

du> A = -Y"b/\" B -w^aw 3 . (2.13) 

B 

These equations allow us to define the connection 1-forms u)\ and u> B if the forms 
duj 3 and du) A are known. Since the 1-forms u> 3 and ui A are the basis forms, then 



(i,j = 1,2,3). (2.14) 



Comparing the equations l|2.8|l and (|2.9|l with the equations (|2.12|) and l|2.13|l . we 

obtain 

u, 3 = -u A = e-^* A u*-e-^ii A ^ A + \Y^~ >lA ~' iB QABU B , (2.15) 



2 

B 



lj a = -u B = ~-e^^-^Q AB u; 3 + e~» B fi A:B u> A - e -^^ B:A u: B , (2.16) 



1 
2 

where 



Qab = <L4:B - Qb-.a, (2.17) 
*A = ^-.A + qA,3- (2.18) 
From (|2.15|) and 12.16fl for the different connection forms we have 

u\ = e-V.sW 1 - \e*»-^-^Q X2 u 2 - e^^w 3 , 

= -ie^-^-^Qa^^e-^^s^-e-^^a; 3 , (2.19) 

^2 = e-^mW -e-^(i 2 W -^-^-^Q 12 u> 3 . 
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Further, in order to culculate the components of the Riemann tensor from the 
second Cartan structure equation 



f 



~R) kl u k t\u l = 0} = du) + wjAw), 



(2.20) 



it is necessary at first to calculate the exterior derivatives of the connection forms 

ens). 

Lemma (Chandrasckhar 8 ). If F is an arbitrary functions of the arguments x 1 , x 2 
and x 3 , then 

d(Fu: 3 ) = ^e-^-» A V A {Fe^)u A ^ + )-^Fe 4 '-» A -^ B Q AB u A i\u) B , (2.21) 

A A,B 



d(Fu) A ) = e-^ A -" B (e» A F) :B u B hu> A + e^"^ (e w f), 3 u 3 Au A , (2.22) 

B 

where Da is an operator, the action of which on an arbitrary function f(x , x 2 , x 3 ) 
is defined by the following expression 



V A f = f.A + q A , 3 f = f,A + (<Za/),3- 
Using this lemma, we obtain 



(2.23) 



du\ = -Y,^'-^A(le 2 ^^-^Q 12 )u A AL 



-l/>—l*l(pl*l—l*2 



(e" 2 /ii :2 ),3U>W + e^>-^(e^-^V2:i),3U 2 Au> 



.2 . , .3 



- U3 [\<jj 



e-K-"* [(e" 2 /i 1:2 ) :2 + (e" 2 -"V 2:1 ) :1 ] L (2.24) 



e~ 2, Vi,3M2,3 - -e 2 ^ 2 ^- 2 ^Q 12 Q 



21 



1 ? 
u> aw - 



1 ^ 

w A" - 



+ 



e «-«-3WQ 12$2 + e-^-^*l/X 2 ,3 



u^Aw, (2.25) 



(e^^ 2 ,3):i + (2e^ 2 Q2i):2 |> , (2.26) 
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Ql2Ml:2 



LO 1 l\tjj 2 - 



+ 



u^A^ 3 - 



w 2 au 3 , (2.27) 



dwj = 



A 

e~ 



(e^-^ m> 3): 2 + (-e^Qi2):i 



, (2.28) 



Tp-2fl! -fj.2 



H2:lQ 



21 



1 2 



e~ 2fl >i:2*2 - -e 2,/ '- 2 ^- 2 ^Qi2Q 2 i 



^2:1* 



2^^ 



u^Ao; 3 - 



2 M2,3<2l 



w 2 Aw 3 . 



Further, from the equations (|2.14l) and l|2.20|l we obtain 



-R\ kl u> k l\J = n\ = &u>\ - wJaw 



(2.29) 



(2.30) 



Substituting H2.24(l - (|2.25(l into this equation and collecting the coefficients at u) k f\u) 1 , 
we obtain the components R\ k i of the curvature tensor. For example, with the object 
to calculate the component i? 212 we must collect the coefficients at w^w 2 in the 
expression for Analogously, the components -R 213 ano - -^223 are obtained from 
Vl\ via t ne comparison of the coefficients at u'au 3 and u 2 Au 3 . In like manner 
from the equation 



1 



Rl kl u k r\J =^1= da; 2 - u\su)\ 



(2.31) 



and equations (|2.26|l - (|2.27Jl we obtain the components i?§ 2 3 an d -R313- Analogously, 
from equation 



1 



(2.32) 



and equations (|2. 28(1 - (|2. 29(1 we have the component i? 313 . Finally, we have the 
following six essential components of the curvature tensor: 



Wl2 e 



^ '' , "^ 2 /il: 2 ):2 + (e W - W M2:l) : 

1 



e" 2, Vi,3M2,3 + ^- 2 ^~ 2 ^Q 12 Q 2l , (2.33) 



i4 3 = -e^-^ 1 {\/2e 2 ^-^-^Q 12 ) - e^-^ (e^>i:2),3+ 

+ e^"^/^ + ie^-^-^^iQai, (2.34) 
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#223 = -e~^-^'D 2 {\/2e 2 *-^-^Q 12 ) - e^-^e^-^iXs- 



ie*-" 1 - 2 " a Qi2*2 - e-^"" 1 *!^, (2.35) 



i?2 23 = - e -*-ftDj( e ^fj) - e-^-^(e^-^/i 2 , 3 ),3- 



e -2 " 1 *!/^:! + ie 2,/ '- 2 ^ 1 - 2 ' 12 g2 2; (2.36) 



i?2 13 = -e-^D^^ftfj) + e -^-Mi(i/ 2e ^-M2Q 21 ) 3 _ 

- e -fi-» [l/2Q 12/ i 1)3 - ^x/ix-a] , (2.37) 



i4 3 = -e-^-^DiCe* - " 1 *!) - e^-^e^'V.sXs- 

- e-^/iiatfa - Je^-^-^QisQai- (2.38) 

2.2 Solutions of equations of geodesic deviation in 
the three-dimensional space 

Let us consider a particular case (/ii = <7i = 0) of the metrics (|2.1|) . In this case 
the metrics (|2.1f) is coincide with the three-dimensional metrics considered in [Sj if 
suppose 

A(x, y,z) = - (e 2 ^ + q 2 e 2 ^) , B(x, y, z) = q 2 e 2 ^ , 

C(x,y,z) = -e 2 ^. (2.39) 

At the condition [i\ = qi = the covariant and contravariant vectors 12. 211 - 12. 31) 
take the form 

e (1)i = (0, 0, -1), e (2)i = (0, -e" 2 , 0), 

e ( a)i = (-e^, 0, g 2 e^); (2.40) 
efo = ( 92 , 0, 1), e l (2) = (0, e"" 3 , 0), 

e| 3) = (e"* 0, 0). (2.41) 

It is easy to see that in this orthonormal basis for the components of the curvature 
tensor we have 

Rjkl= -Rijkl- (2-42) 

It is well-known that the Riemann tensor Rijki has the following symmetry 
properties: 



(2.43) 



It is easy to show that the symmetry properties H2.43|l decrease the num- 
ber of independent (essential) components of the Riemann tensor from n 4 to 
n 2 (n 2 — 1)/12, where n is a dimensionality of the space. In the case of three- 
dimensional space we have six independent components of the curvature tensor: 
#1212, #1213, #1223, #1313, #2313, #2323- Further, using H2.42j)-I|2.43j). we see that 



Rijkl 


= Rklij, 


Rijkl 


= —Rjikl 


Rijkl 


= —Rijlk 
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in the system ()1.4f) among the four components of the curvature tensor only three 
are independent, namely, i? 2 2 i, -R131 and i? 131 (or i?i 2 i)- The latter two components 
are coincide with each other in virtue of (|2.42[l - il2.43|l . Therefore, 

p3 p2 

rt 121 — - fl 113- 

Hence it immediately follows that the conditions (|1.8(l and the system l|1.7|) are 
reduced to the form 

A 2 - R 2 = -u 2 , X 2 - R* = -u 2 , 
121 131 > 

R'l21 = — -<H13 = U x'i 



1 0\ d 



11 dx 2 \-u x -u 



u , 



,2 



It is easy to see that the matrix equation l|2.45[l corresponds to the Zakharov-Shabat 
system [5] 

— \-lXlpl = U1p 2 , 

OX 

^p--i\ip2 = -Ufa. (2.46) 
ox 

Thus, in the orthonormal basis \2.40^ - \Z.41}) , related with the metrics at the 

condition = q\ = the AKNS-system for the equations of geodesic deviation 
is reduced to the Zakharov-Shabat system. Moreover, instead the two potentials in 
AKNS-system we have now only one potential in ZS-system. 

Let us calculate the independent components of the curvature tensor in the 
system (|1.4|) for the metrics (I2.1|) at the condition fj, x =q r = 0. From (|2.33|) . (|2.34l) 
and iflHSj) we have 

^212 = -^ii-zxlx-je-^^^e^ + l), (2.47) 
^13 = -e^-^.ii + |e^V,i92,i - ie*-*«/i2,i«a,i. ( 2 - 48 ) 
^313 = -^11-^1 + \e 2 ^ 2 ^ q 2 A . (2.49) 

So, in the case of the metrics (I2.1|) our problem of solving of the equations of 
geodesic deviation is reduced to the Zakharov-Shabat problem (|2.46f) . It is known 
that fundamental solutions (Jost functions) of ZS-problem are defined by the fol- 
lowing expressions 



Vi{x,\) = e- lXx + J dx'Aiix,^- 1 ^' , 

— CO 

X 

<fi2(x,X) = { dx'A 2 {x,x')e~ lXx '; 



(2.50) 



ift(x,\) = J dx'Bx(x,x')e iXx \ 

X 

CO 

(pt(x,X) = e lXx + l dx'B 2 {x,x')e lXx . 



(2.51) 
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These solutions are linearly dependent: 



if (a, A) = cn(X)ip + (x, A) + ci 2 (\)<p + (x, A), 
<^ + (x,A) = c 2 i(\)(p~(x,\) + C 22 (\)lfi~(x,\), 



(2.52) 
(2.53) 



where 



M-Hfz*)- ^>-(*sr-V (2 - 54) 

Further, the pair of Gel'fand-Levitan-Marchenko integral equations can be derived 
from (|2.52() by means of the Fourier transform: 



A 2 (x,y) + fl L (x + y) + / dx 1 A 1 {x,x')9, L {x' + y) = 0, 



A 1 (x,y)+ / dx'A 2 (x,x')n L (x' + y) = 0, 



x > y. 



Analogously, from l|2.53[) we have the following pair 

oo 

B 2 (x,y) + { dx'B 1 {x,x')^ R {x l + y) = 0, 



(2.55) 



where 



-Bx{x, y) + n R (x + y) + J dx'B 2 (x, x')9. R {x' +y)=Q, 

X 

x < y, 



N 



Ml = r L (z)-zV" T 
* — * i ■ 



(=1 

N 



C22(A;) ^ 

cia(Aj) 



-iX t z 



(2.56) 



(2.57) 
(2.58) 



^c 2 i(A z ) 

Thus, the potential u(x) is expressed via the kernels A\, A 2 and B\, B 2 as follows 



u = —2A 2 (x, x), u = —2B\(x. x), 
2 _ n dAi{x,x) 2 _ n dB 2 {x,x) 

LL — ^ 7 , Li — ^ 7 

ax ax 



(2.59) 



In the case of a reflectionless potential (tl(z) — 0) the system of Gel'fand-Levitan- 
Marchenko integral equations may be solved explicitly. In this case we obtain 



N 



n \ - c 22 {\p _ 



■ikiz 



(2.60) 



The system <|2.55(1 is reduced to algebraic equations and the potential u(x) is defined 
by a following expression 



u = 2— arctan 
dx 



Imdet(7 - iM) 



Redet(Z-iM) 



(2.61) 
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where 

im Lj i(lii+K . )x 



MlJ = J , 1 e-* Ki+K *> x , (2.62) 
m Lj = 



_.C 22 (\j) 



ci 2 (Aj)' 

Here «j are the poles of the transmission coefficient Xj,(A) = - . Analogous 

relations take place in the case of system <|2.56[1 . 

Let us return to the equations of geodesic deviation. From the conditions on the 
curvature tensor (|2.44ll it follows that 

p2 p3 

rt 121 — - tl 131) 

^121 = u x- 

Substituting the expressions 12. 4711 - 12. 49fl into the latter equations, we obtain 

V',ii + ^ 2 i-/i2ai-M2,i-^^ 2 9 2 2 i(e 2V ' + ^) = 0, (2.63) 
\e^'^q 2 ,ii + ^e^ 2 VAi92,i - ^-" 2 ai 2 ,i«2,i = u,i. (2.64) 

Thus, we have the system of differential equations (|2.63|) - (|2.64|) as the conditions 
on the potential u{x). The explicit form of u{x) we will find by means of the 
inverse scattering problem. Moreover, the potential u(x) depends parametrically on 
variables y and z. According to widely accepted methods [6l ll3llTl] . the dependence 
on variables y and z may be represented by a nonlinear integrable equation. Indeed, 
the dependence on y for ?pi and t/> 2 from I2.46fl may be expressed in general form 

ipi y = Aipi + Bifci, 

= CiPi+D^. (2.65) 

The compatibility conditions of (|2.65() with i|2.46[l give (at this point, X y = 0): 

A x = u(B + C), (2.66) 
B x + 2i\B = -2uA + u y , (2.67) 
C x -2iXC = -2uA-u y , (2.68) 

here D x — —A x . Further, let us suppose that A = ^jj a„A", B — Y^o^n^™ 
and C — ^qC„A™. Substituting these series into 12.67|I - I|2.68J) . we obtain for the 
coefficients a n , b n and c n the following expressions 

a.3=a 3 (y), b 3 = c 3 = 0, 
a 2 =a 2 (y), b 2 = -C2=ia 3 u, 

1 2 . 1 1 

fl i = ~2 fl3M ' 1 = ~2 a3Ux + m2U ' ° = ~2 a3Ux ~ m2U ' 

1 i 1 

a = --a 2 u 2 , b = c = -a 3 (u xx + 2u 3 ) - -a 2 u x . (2.69) 

In the equations (|2.67l) - l|2.68[) the components independing on A give the evolution 
equation u y = b$ x + 2aou. Using the obtained above expressions (|2.69ll . we obtain 
for the coefficients ciq and b : 

u y = -jia 3 (6u 2 u x + u xxx ) - a 2 (u 3 + ~u xx ). 
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Suppose ai — and 03 = Ai we have the modified Korteweg-de Vries equation 

u y + 6u 2 u x + u xxx = 0. (2.70) 

Thus, the dependence on parameter y for the potential u is defined by the mKdV 
equation. Thus, the system 12.65( 1 has a form 

iply = 2i\(u 2 - 2A 2 )f/'i + (4A 2 u + 2i\u x - 2u 3 - u xx )tp 2 , 

ip 2y = (-4A 2 m + 2i\u x + 2u 3 + U xx )ipi - 2i\(u 2 - 2X 2 )-ip 2 . (2.71) 

Solutions of the modified Korteweg-de Vries equation can be found by the standard 
procedure El E] ■ When u->0we see that the dependence on y is described by 
a limiting form of the equations (|2.71|l 

iply = -4iA 3 ^i, 

^ 2y = 4i\ 3 ^ 2 . (2.72) 
Let us assume that ip — (^jf^ * s proportional to the fundamental solution (p~ 

at x — > —00. Then, at x — > —00 we have ip{x,y) = f{y)<p~ — > }{y)^ lXx i^j ■ 

Substituting ipi = f{y)e~' lXx into the first equation from (|2.72|l . we obtain f(y) = 
/(0) exp(— 4iA 3 y). From 12.52|l at x — > +00 it follows that 



V> = f(y)<p- — /(o) e - 4 ^ 



Cll (A,y)e^ Q+c 12 (A,2/)e- aa: Q 



(2.73) 



Substituting it again into (|2.72|l . we obtain that cyi = 0, en = 8iA 3 , whence 

cu[\,y) = ci 2 (A,0), 

Cll (X,y) = Cll (A,0)e 8 ^. (2.74) 

The analogous calculations for ip ~ ip + give 

C22(A,2/)=c 22 (A,0)e- 8lA ^. (2.75) 

Using the dependence on parameter y given by the relations 12.74fl - (|2.75l) , we have 
for 12.62(1 and (|2.57l) the following expressions 

m Ll (K U y) = -t C22i . Kh l\ = m Ll ( Kl ,0)e- StX >y, (2.76) 

Cl2(«J,0) 

OO 

Vl(z, y)= f ^Rl(K 0)e- six3 «- a * + ]T m Ll ( Kl ,0)e^y~^ z , (2.77) 

where i?z,(A, 0) = — c 22 (A)/c 2 i(A). Further, from (|2.59|) it follows that the potential 
u(x,y) is expressed by the kernel of Gel'fand-Levitan-Marchenko equations (|2.55|l 
as 

u{x,y) = -2A 2 (x,x,y). 

In case of the reflectionless potential (Rl(\, 0) = 0) the integral equations (|2.55|) 
are solved explicitly. In this case, the potential u(x,y) is defined by the formula 
((2.61(1 with the matrix M of the following form 

My = j m ^( /8 J'W) e -<(m+Hi)x < (2 78) 

Kj + Kj) 
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In the simplest case of the one-soliton solution (N =1), the matrix M is reduced 
to a scalar M = «(toi/2ki) exp(— 2inix). Taking into account the relation l|2.7tj[l . 
we see that the matrix M at k = iX can be written as 

= mi(0) 2AB _ 8A 3 M 
2A 

Thus, in case of the one-soliton solution, the potential, defined by the expression 
(|2.61(l with the matrix (|2.78() . is reduced to the form 



/ \ „ d 
u(x, y) = —2— arctan 
ox 



2A 



(2.79) 



or 

u(x, y) = ±2A sech (2Aa; - 8A 3 + 5), (2.80) 

where S = In [mi(0)/2A]. For mi(0) < we take the upper sign, for toi(0) > the 
lower sign. 

Supposing the analogous dependence on the parameter z, that is, defining it by 
the modified Korteweg-de Vries equation of the form 

u z + 6u 2 u x + u xxx = 0, (2.81) 

we came in the case of the one-soliton solution to the following dependence 

u(x, y, z) = ±2Ascch(2Aa; - 8A 3 y - 8A 3 z + 5). (2.82) 

Thus, we see that dependence of the potential u on the parameters y and z is given 
by the mKdV equations (|2.7l)|l and (|2.81|l . Let us consider now how the vector 
of geodesic deviation rj may be expressed via the fundamental solutions (p^ of the 
Zakharov-Shabat problem (|2.46|l . We will consider here two particular cases of the 
system lj2~fi31l - lj2~Mll . 

2.2.1 V = 

In this case, the coefficients (I2.39|) of the metrics (|1.3|1 are 

A = - (e 2 ^ + q 2 2 ) , B = q 2 , C = -l. (2.83) 
And the system (|2.63|) - (|2.64l) is reduced to the form 

/"2,n + M2,i + l^ll,! = 0, (2.84) 
-e-^q^u ~ ^e-">2,i92,i = u, x . (2.85) 
The latter equation, obviously, may be written as 

whence 

u=\e-^q %l . (2.86) 



2 

On the other hand, in virtue of H2.82|l we have 
1 



q 2A = ±2Asech(2Aa; - 8A 3 y - 8A 3 z + 5). (2.87) 
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Using the well-known relation sinh 2 A — 2 cosh A sinh A, we can write (choosing the 
upper sign) l|2.87fl in the form 

= 8Xsmh(2Xx-8X 3 y-8X 3 z + 6) 
924 sinh(4A.T - 16X 3 y - 16A 3 z + 26) ' 1 ' ' 

Whence, supposing 



sinh(4Aa; - 16A 3 y - 16A 3 z + 25) ' 
g 2 ,i = 8Asinh(2Ax-8A 3 y™8A 3 z + <5), 



we obtain 



or 



^ 2 = lnsinh(4Ax- 16A 3 y - 16A 3 z + 2<5), (2.89) 
q 2 = 4Acosh(2Ax-8A 3 y-8A 3 z + ^), (2.90) 

^ 2 = msinh(4Ax- 16A 3 y - 16A 3 z + 2<5), (2.91) 
q 2 = -4A cosh(2Ax - 8A 3 y - 8A 2 + 5) (2.92) 

in the case mi(0) > 0. 

Thus, using H2.50fH2.51fl an d l|2.59l) . we obtain that solutions of the equations 
of geodesic deviation (|1.4fi in the case of the metrics <|2.83fl are expressed via the 
fundamental solutions of the Zakharov-Shabat problem as follows 

X 

rj 2 ~ <Pi(x,X) = e~ lXx + i j dx 1 j u 2 dxe~ tXx ' , (2.93) 
t? 3 ~ <to(x,\) = ~ j dx'ue~ iXx ' rim) 



~ (ft{x,X) = -i / dx'ue lXx \ (2.!). r ,i 



rf ~ vt(x,\) = e tXx - i / dx' / u 2 dxe iXx , (2.!)<il 



where 



u=-e^ 2 <Z 2 ,i, (2-97) 



and the functions /i2 and 92 are related by the equation 

M2,ii+M2,i + \e- 2 ^ql 1 = 0. (2.98) 

For example, in the case of the parameter dependence on y and z described by 
the mKdV equations (|2.70fl and l|2.81fl at N = 1 (onc-soliton solution) we have the 
following integral representations 

X 

rj 2 ~ e~ lXx + X J dxtanh(2Xx~8X 3 y~8X 3 z + S)e- lXx , 

— OO 

X 

rf ~ =fA J dxsec\i{2Xx-8X i y-8X i z + 8)e- iXx (2.99) 
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r, 2 



TA J dxsech{2\x - SX 3 y - 8A 3 z + S)e tXx , (2.100) 

X 

OO 

„ e lXx - X [ dxta,nh{2Xx-8X 3 y-8X 3 z + S)e lXx 7 (2.101) 



where for mi(0) < we take the upper sign and the lower sign for mi(0) > 0. The 
constraint l|2.98|l gives (here the functions fi 2 and q 2 are defined by (|2.89|) - (|2.90ll or 

ESD-ES3): 

cosh 2 (4Ax - 16A 3 t/ - 16A 3 z + 25) + 3sinh 2 (2Air. - &X 3 y - 8A 3 z + 5) = 1. 
2.2.2 ^2 = 

In this case, the coefficients (I2.39|) of the metrics (|1.3ll are 

A(x,y,z) = -(l + B(x,y, z) = q 2 e^ , C(x, y, z) = -e™. (2.102) 

The system (|2.63[) - (|2.64[) is reduced to the form 

^ 11 +^l-lql x (eW+ 1 -) = Q, (2.103) 

ie^^n + |e^,i«2,i = (2-104) 

lb 

Using the substitution 9 = -Sj, we obtain from the latter equation 

~ (e 9/ W) x = u,i. (2.105) 

Therefore, in this case the vector of geodesic deviation r] is also expressed via the 
fundamental solutions ip T (x, A) of the form (|233>(|2~§4"|) or ((235^) - ((23^ . At this 
point, 

u = ^e 9 / 3 (7 2il (2.106) 
and the functions 9 and 52 are related by the equation 

f,n + ^i-|ft,i(e a/!W + i)=0. (2.107) 

In the case of the one-soliton solution, a potential u is defined by the expression 
(|2.82|) , and the functions 9 and 92 are respectively equal to 

6 = 31ncsch(4AiE- 16A 3 y - 16A 3 z + 2(5), (2.108) 
q 2 = ±cosh(2Aa;-8A 3 ?/-8A 3 z + (5). (2.109) 

More complicated case \i 2 7^ 0, ^ ^ and also multi-soliton solutions will be 
considered in a separate paper. 

3 Four-dimensional space 

3.1 3x3 matrix Schrodinger equation 

Let us consider the following linear problem 

ip A = iCDip + Nib, (3.1) 
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where ip 4 = 77^- "0? x — it'i C is a spectral parameter and -0 is a 3 x 1 matrix 
ox 

(vector) of the form 







f 


#12 


#13 




#21 





#23 




l#31 


#32 






The 3x3 matrices D and N (a potential matrix) are 



D = 



The system (|3.1f) may be rewritten (see Appendix) in the following form (3x3 
matrix Schrodinger equation) 

-Iip,44 + V\ip = d 2 C 2 ^, (3.2) 

where / is a 3 x 3 unit matrix, d = (d\, d2, d 3 ) and 

/# 2 #>i + N 13 N 31 N 12A + N 13 N 32 N 13A + N 12 N 23 
91= \ N 21A + N 23 N 31 N 21 N 12 + N 23 N 32 N 23A + N 21 N 13 
\ #31,4 4- #32 #21 N 32A + N 31 N 12 N 31 N 13 + N 32 N 23/ 

Further, it is easy to see that for the metrics 1)1. 9[) the geodesic deviation equation 

frmji . 

Vm + RIuV 1 + RI24V 2 + RI34V 3 = 0, 

V% + RluV 1 + R 2 42aV 2 + RI34V 3 = 0, (3.3) 

v%4 + Rluv 1 + R424V 2 + RI34V 3 = 0, 

can be rewritten in the form of the 3x3 matrix Schrodinger operator 

-VM + i-Rlu + dieW + i-R^f + i-R^v 3 = d?cV, 

-V 2 44 + (-R 2 414W + (-R 2 242 + d 2 2 ( 2 W + (-Rl 3 4)V 3 = ^CV, (3-4) 
-V%4 + {-RlliW + (-R 3 424)V 2 + (-R 3 434 + dfcV = 4CW ■ 

Comparing these equations with equations (|3.2|l , we obtain the following conditions 
on the curvature tensor: 



die 


R444 — 


N 12 N 2 i 4 


" #13#31 


dlC 2 


~ R 2 24 = 


#21 #12 4 


- #23#32 


die 2 


^434 — 


N 31 N 13 4 


- #32 #23 


— R-424 


= #12,4 4 


"#13 #32, 




pi 

rl 434 


= #13,4 4 


- #12#23j 




— -R414 


= #21,4 4 


-#23#3i, 




E>2 

— -"-434 


= #23,4 4 


-#2i#i3, 




-n.414 


= #31,4 H 


-#32 #21, 




R 3 

— -n.414 


= #31,4 H 


-#32#21- 





(3.5) 



3.2 Chandrasekhar metrics 



In the 4-dimensional space with a signature (—,—,—,—) the Chandrasekhar metrics 
is defined by the following expression [5J 

ds 2 = -J2 e2tlA (dx A f ~ ^(dx 1 - <lAdx A ) 2 , (3.6) 

A A 
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where A = 2,3,4. ip : [ia and qA are the functions on variables x 1 , x 2 , x 3 , x 4 . 

The orthonormal tetrad, related with the metrics (|3.6|l , is defined by the follow- 
ing covariant basis vectors: 



3(4)i 



(-e^, 0, 0, 0), e (1)j = (g 4 e^, -e* g 2 e^, q 3 e^), 



e(2)» = (0, 0, 



0), 



j (3)i = (0, 0, 0, -e" 3 ). 



(3.7) 



And also the contravariant basis vectors are 



'(4) 

,i 

'(2) 



(e^ 4 , <? 4 e~'' 4 



(0, g 2 e-^, e"" 2 , 0), 



0, 0), = (0, e-*, 0, 0), 



'(3) 



(0, fie""', 0, e-«). 



From (|3.7|l and (|3.8|) it is easy to see that 



(3.8) 



e {a) e {b)i = V(a)(b) 



-10 

0-100 
0-10 
-1 



Therefore, in this orthonormal basis for the components of the curvature tensor we 
have always 



D» 



-Ri 



hi I: 



(3.9) 

Moreover, at dt = — idx 4 , v = /X4 and cj = ^4 there exists an analytic contin- 
uation of the basis l|3.7|) - (|3.8|) with the signature (— , — , — , — ) onto a basis with a 
signature (— , — , — , +), the covariant and contravariant vectors of which have the 
form 



e(i)i 
e (3 )i 
e (i) 

e (3) 



we' 



g 2 e^, g 3 e^), e (2)i = (0, 0, 



-,/<2 



0), 



= (0, 0, 0, -e^), 
(0, c-* 0, 0), 
(0, (fee"" 3 , 0, e-«), 



e (4)i = (e", 0, 0, 0). 



'(2) 
'(4) 



(0, g2 e^ 2 , c-"», 0), 
(e"*, we'", 0, 0). 



(3.10) 



It is obvious that in the orthonormal basis H3.7[) - l|3.8|) among the nine compo- 
nents of the curvature tensor of the system (|3.3|l only six are independent, namely, 

pi p2 p3 pi pi p2 
-"■414; /l 424; ^434: /l 424) ^434' ^434 • 

In the orthonormal basis H3.7fl - H3.8fl f° r the metrics l|3.ti[) these components have the 
form 



- fli4i4 = -e'^-^Di (e^ 4 ^) - 

- e-*"*** (e^-^ M4 ,i) x + ~ 



^ 2 ^ 4:2 



-2^3 



*3M4:: 



,2V-2,U4 r^-2M2Q2 



24 



3 -2M3n 2 



^34 J ) 



(3.11) 



i? 2 



e 2/i3 /i4:3M2:3 ^ 



"M4:2 



; 2 4 



-2(/' 



M2,iM4,i) (3-12) 



3434 



-M3-M4 [(gW-W^) + ( e M4-M3 M4:3 ) 3 



6 2/i2 /i4:2M3:2 - ^6 



2ip-2fi 3 —2fj, 4 ^2 
^34 



e- 2V, ^ 3 ,i^4,i, (3-13) 
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#1424 = e^ 2 ^-^Q 42 (* 4 - ^2:4) + X -e-^-^ (e^-^Q 42 ). A + 

+ ^- 2 ^^Q 32f i 4:3 + e~^-^ (e-^V4,i) :2 - e-^> 2 ,iM4:2, (3.14) 



#1434 — e 



1 



+ ^e^- 2 ^-^g 23 ^4:2 + e-^-^ (e-^ M 4,i) ;3 - e-^ 3 M 3 ,iM4:3, (3.15) 



#2434 = e A ' 2 M3 [^4:32 + M4:3 (^4 ~ M3) :2 ~ ^4:2^2:3] ~ 

- \e 2 ^^~ 2 ^-^Q Z4 Q 42 - ie-^-^Q 32 /i 4 ,1, (3.16) 

where 

2>a/ = f,A + {q A f)^, 

= 1p:A+qA,l, 
QAB — q A -.B — QB-.Ai 

f:A = f, A +q A f,i- 

3.3 Solutions of equations of geodesic deviation in 
the four-dimensional space 

It is easy to see that the Chandrasekhar metrics (|3.6H coincides with the metrics 
(|1.9|) at /14 = q 4 = 0. In this case, the orthonormal basis is reduced to the form 

e (1)< = (0, -e* q 2 e^, g 3 e^), e (2)i = (0, 0, -e» 0), 

e (3)i - (0, 0, 0, -e»), e (4)i = (-1, 0, 0, 0); 

e* (1) = (0, e-*, 0, e-^), e * (2) = (0, q 2 e^\ e~^, 0), 

e* (3) = (0, gae-"», 0, e"^), e l (4) = (1, 0, 0, 0). 

It is obvious that in this basis we have i?™ fei = —Rijki, and the components of the 
curvature tensor (|3.11JI - I|3.16[1 are 



-#2424 = -M2,44 - /i2,4 - Te 20 ~ 2M2 g 2 ,4, (3.18) 



#1414 = -V,44-^f 4 + \^ [e^ 2 92,4 + e- 2 ^ql 4 ] , (3.17) 

3 
— < 

4 

-i?3434 = -M3,44 - /il 4 - J^-^g 2 4j ( 319 ) 

#1424 = - 2 e^"" 2 [92,44-^2,432,4 + 3^,452,4], (3.20) 

#1434 = -^-^ \q ZM - M3 , 4 g M + 3^ 49m ] , (3.21) 

#2434 = ^e 2 ^-' l2 -^g 2 ,493,4. (3.22) 

Further, let us define now the evolution equations related with the problem (|3.1|l . 
We consider the following system 

iP A = i(Dip + Nip, 

tp,i = Qi>, 
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where w % = ttt, x is a parameter of the considered problem. D, N, Q be the 3x3 
matrices. At this point, D is diagonal: D = diSij, di — const; N is such a matrix 
that Na = 0. From the compatibility condition ip A 4 — ij) A i and the requirement 
C.i = we obtain 

Q A = N A +iC[D,Q] + [N, Q\. 
Decomposing Q in the form 

Q = Q (1) C + Q (0) , 

we have = N i + [N, Q^], whence we obtain the system of nln — 1) equations 
(see©): 

Nijs - aijNi jA = y^(q;fc - a k j)Ni k N kj , (3.23) 

k 

where 

1 qi - qj 

i di — dj 

Equations (|3.23[) may be reduced to the standard system of nonlinear equations of 
three-wave interaction. Namely, we obtain 

<9i,i + CiQ M = ijiQiQt, 

Q2,i + C 2 Q 2A = h 2 Q*iQ* 3 , (3.24) 
Q3,i + C 3 Q 3A = hzQ* x Ql, 

where 717273 = —1, 7i = ±1 and 

N 12 = -iQs/VNsR, N 31 = -iQ 2 /VW^, 

N 23 = +iQi/VPl 2 ~Pw, N 13 = -7173^31= ( 3 - 25 ) 

N 32 = 7273 ^23' N 2l = 7172^1*2, 

here 

.C\C 2 C 3 

qj = -i — -p-, , Pij =di- dj = Cj - Ci, 

C 3 > c 2 > c\. 

In the system Ij3.24|l there is a decay instability (for the waves with positively defined 
energy) if the sign of one j n is different from the other, and also there is an explosive 
instability when 71 =72 = 73 = —1- Solutions of the system l|3.24|l was obtained 
by Zakharov and Manakov in 1973 |10HlllH2*] . They have the form 



e Xi(^ 4 -C , 3 a; 1 -¥'3_ 



7172 C ' "^ e-* 1 ^ 4 -^ 1 -^) 

d* - C3 



(3.26) 



-4XlX3/3l37273 - l ^-C 3 x 1 ^ 3 )- l ^(x i -C 1 x 1 -^ 1 ) /o 97 ^ 



Q 3 = VPAlil2 2 -^-e l ^~ c ^~^ 



7273 £i Ge-Xai^-Cix 1 -^) 



Cl*-C3 



, (3.28) 
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where 



S 



e Xl(» 4 -C3X 1 -y3) _ ^, l ^ 2e -Xi(^ 4 -C32: 1 -<p3) 



7 73 (c"i-a)(cr-c 3 ) ' ( j 

4i — — 3 , S3 ^ 



/J 



2:1 



Supposing now that the matrix AT is real (TV* = N) and choosing 71 = 73 
1, 72 = —1, we obtain from (|3.25|1 



N 12 = -Re(iQ 3 /VPl^), N 31 = -Re(iQ 2 /Vft^3~), 

N 23 = +Rc(iQ 1 /^ 2 ~^), N 13 = -N 31 , (3-30) 

^32 = -N 23 , N 21 = -JV12. 

It is obvious that the latter three conditions in H3.30|) are equivalent to antisymmetry 
of the matrix N. 

Thus, we assume that the potential matrix N is real and antisymmetric. Taking 
it into account and also the expressions (|3.17|) - (|3.22|) . we obtain from the conditions 
on the curvature tensor (|3.5() the following system of differential equations: 

-^,44 - ^4 + [e- 2 ^ql A + e" 2 ^] = ^2 + *?3 + d\? , (3.31) 



- M 3,44 - M!,4 ~ I^^qIa = ^13 + ^23 + 



^ M2 [<Z2,44 - M2,4<?2,4 + 3^,452,4] = N 23 N 31 , 

>-Ms [g 3 44 _ fJj3A q 3 4 + 3ip A q 3A ] = ^12-^23, 



,2i/'-M2-A'; 



<?2,4<73,4 = N 31 N 12 . 



(3.32) 
(3.33) 
(3.34) 
(3.35) 
(3.36) 



Obviously, this system has a great number of particular cases. For example, let us 
consider one simplest case. 

3.3.1 q 3 = const, ip = [i 2 = 

In this case for the metrics l|1.9[l we have 



3ii = -1, 322 = -(1 + q 2 ), 933 = ~{e 2f23 + const 2 ), 
.912 = 2g 2 , 3i3 = 2const, 323 = -2constg 2 - 
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The system (|3.31[) - (|3.36[) is reduced to the form 



ll 4 = #12 + #13 + 


die 2 , 


(3.37) 


|?2 2 ,4 = ^l 2 2+^2 2 3 


+ d 2 { 2 , 


(3.38) 


M3.44 - A*3,4 = #13 




(3.39) 


ig 2 ,44 = #23#3i, 




(3.40) 


N 12 N 23 = 0, 




(3.41) 


N 31 N 12 = 0. 




(3.42) 



From the latter two equations it follows that Ni 2 = 0. Therefore, in this case, the 
potential matrix N has a form 

-7V 3 i N 
iV 23 
-iV 23 

Further, at d\ = -4/3d| from IpOTfl - lpO^fl it follows that 



"2 92,44 = #23#31, 



"M3,44 - ^3,4 



-g^ + ^C 2 , 



#13 + = 0. 



In accordance with l|3.30|l . the components iV 23 and N 3 i are defined as 



N: 



23 



2X3 
D 



2(^23X1 - ^12X3) 



x cos 6 (a; 4 - C lX x - (p x ) e -xx(**-c 9 « 1 - Va )_ 
— 2 sin £3 (cc 4 — Cix 1 — ifi) coshxi (x 4 — C3X 1 — ^3) 



(3.43) 
(3.44) 
(3.45) 



(3.46) 



#si = - 



4XlX3/5l3(/323^1 - 



(foatl ~ /3l26) 2 + (/?23Xl " /?12X3) 2 £> 

sin^i (x 4 - C3X 1 - (p 3 ) cos£ 3 (a; 4 - Cix 1 - (pi) + 

+ cos£i (a; 4 - C3X 1 - ip 3 ) sin £3 (x 4 - C\x x - (pi) 



, (3.47) 



where 

D = 4coshxi (x 4 — C 3 x x — (p 3 ) coshx 3 (a; 4 — Cix 1 — ipi) + 

, 4 XlX3 

PlM + Xl) + 2/3 12 /3 2 3(^3 + X1X3) + (3 2 2 M + Xl) 

After very cumbersome but elementary calculations it is easy to verify that solutions 
of the system H3.43M3.45H exist. 
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Appendix 



Let us consider the system (|3.1|) with the matrix 



f-di 
D = d 2 
V -d 3 



) 



Differentiating i|3.1[l and excluding the first derivatives ip A , we obtain the following 
system 



- V>i,44 + (N 12 N 21 + NisNsiWi + (N 12A + N 13 N 32 - i{dxN 12 + iQd 2 N 12 )^ 2 + 

+ (N 13A + N 12 N 23 - i(dxN 13 - i(d 3 N 13 )iP 3 = C 2 d?Vi, (A.l) 

- ^2,44 + (iV 21) 4 + ^23^31 + i(d 2 N 21 - i^N^)^ + (N 21 N 12 + N 23 N 32 )<if> 2 + 

+ (N 23A + N 21 N 13 + i(d 2 N 23 - i(d 3 N 23 )Tp 3 = ( 2 d 2 t{j 2 , (A.2) 

- V-3,44 + (iV 3 i,4 + ^32^21 - i(d 3 N 31 - iCdiN 3 i)ipi + (N 32A + N 31 N 12 - iC,d 3 N 32 + 

+ iQd 2 N 32 )ij 2 + {N 31 N 13 + N 32 N 23 )ip 3 = t 2 d 2 3 ^ 3 . (A.3) 

Analogously, for the matrix 



- V-1,44 + {N 12 N 21 + N 13 N 31 )Tpi + (N 12A + N 13 N 32 + iQd x N 12 - iQd 2 N 12 )ijj 2 + 

+ (N 13A + N 12 N 23 + iC,d 1 N 13 + i(d 3 N 13 )^j 3 = fd 2 ^ , (A.4) 

- -02,44 + {N 21A + N 23 N 31 - i(d 2 N 21 + iC,d x N 2X )ij)i + (N 21 N 12 + N 23 N 32 )ip 2 + 

+ {N 23A + N 21 N 13 - i(d 2 N 23 + i(d 3 N 23 )i; 3 = ( 2 d 2 ^ 2 , (A.5) 

- V-3,44 + (^31,4 + ^32^21 + ^4-^31 + K^l N 31 )lpl + (N 32A + N 31 N 12 + l(d 3 N 32 - 

- i(d 2 N 32 )ip 2 + +{N 31 N 13 + N 32 N 23 )ip 3 = C 2 d& 3 . (A.6) 
Adding the systems (|A.ip> - (|A.3fl and (|A.4|I - I|A.6(1 . we obtain in the result 

- Vi,44 + (Ni 2 N 21 + N 13 N 31 )ip! + (N 12A + N 13 N 32 )4> 2 + 

+ {N 13A + N 12 N 23 )^ 3 = ( 2 dl^ , (A.7) 

" 02,44 + {N 21A + ^23^31)^1 + (N 21 N 12 + N 23 N 32 )ip 2 + 

+ (A 23 ,4 + ^21^13)^3 = C 2 ^2, (A.8) 

- 03,44 + (JV 8ll 4 + ^32^21)^1 + (^32,4 + ^31^12)02 + 



It is easy to see that the latter system can be rewritten in the form of 3 x 3 matrix 
Schrodinger equation (|3.2|) . 



ld x 
D = -dk 
\0 d 3 



we have the system 



+ (^31^13 + ^32^23)^3 



( 2 d& 3 . (A.9) 
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